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ABSTRACT 










A simple analytic model is formulated which should prove useful for 
evaluating a satellite surveillance system. This model describes the 
fraction of orbits in which a long lived satellite will "see" a selected 
position on the earth’s surface. The fraction of orbits, denoted by P, 
is a function of the latitude of the ground position in question, the 
orbital inclination of the satellite, and the sweep width of the satel- 
lite’s sensor. 

P is used to find estimates for the mean time to next sighting of a 
selected ground position given that position has suddenly become of in- 
terest, to estimate the number and locations of ground readout stations, 
and to estimate the number of times a satellite would sight a transitting 
object on the earth’s surface. In addition, the application of P to the 
problems of optimum orbit, evasive countermeasures, and satellite weapon 
delivery is discussed. 
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1. Introduction. 



Rapid developments in space technology have resulted in an increased 
use of unmanned earth-orbitting satellites for both civilian and military 
purposes. These satellites provide information on outer space, the earth* s 
atmosphere, and the earth’s surface. Satellites for earth surveillence 
offer a distinct advantage over earth-bound systems because of their abil- 
ity to scan large segments of the earth in a relatively short period of 
time, and because they are unencumbered by national air-space restrictions. 

The Tiros and Nimbus series weather satellites represent well known 
examples of the use of satellites by the United States for earth surveil- 
lence. Tiros III, launched in 1961 was equipped with photographic sensors 
and provided information on storms and storm paths. Nimbus I, launched in 
1965 served a similar purpose but made use of television cameras as sen- 
sors . 

Although weather received the first attention, it is not the only 
phenomenon wDrthy of satellite surveillence. Nuclear test explosions are 
of sufficient interest that two Vela Hotel satellites were launched by 
the United States in 1965. Additional satellites apparently have been 
launched for certain classified military purposes. The Samos series 
satellites evidently have, for the last five years, been photographing 
military installations in Russia and Red China. [ 6j Seven military recon- 
naissance satellites (possibly of the Samos series) were reported to have 
been launched during 1965. M These satellites had polar orbits and 
mission lives of about five days. The observations made by these satel- 
lites are unknown, but could include ship and troop movements or test 
firings of ICBM’s, as well as military installations. 

Given a surveillence task, alternative satellite systems might be 
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designed to perform this task. We may want to decide which system or com- 
bination of systems is most cost-effective. Or we may want to compare 
the alternatives with more conventional systems such as U2 flights. To 
make such comparisons we must be able to evaluate the effectiveness of 
each type of surveillence system. This paper will be concerned with eval- 
uating satellite surveillence systems. 

What elements are important in the evaluation of a satellite surveil- 
lence system? This is a difficult question to answer. The appropriate 
elements must include both the nature of the phenomena to be observed, and 
the available hardware, including the satellite, its components, and ground 
support equipment. The nature of the phenomena to be observed may be such 
that one satellite of a particular type is sufficient. Or, it may be 
found that multiple satellites of one or more types are required. The 
appropriate number will depend upon the effectiveness of each individual 
satellite. Thus, an initial appraisal of a system should be based on an 
individual satellite’s performance. 

What features of an individual satellite’s performance should be ex- 
amined? Two features appear to be dominant. For those phenomena which 
occur randomly in time and location on the surface of the earth, it seems 
desirable to make some statement about the expected proportion of the total 
phenomena of interest the satellite will observe. For other phenomena we 
are more interested in seeing the result of an event rather than the event 
itself. Thus, we would like to know the average elapsed time before a 
satellite observes a location which has suddenly become of interest. (For 
example, if on the average 100 natural disasters occur in a year in the 
area covered by the satellite, then we may want to know the average time 
lapse between occurrence of a disaster and observance of the place of 
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disaster*) 



The present procedure for satellite system evaluation is to construct 
an elaborate model describing the location of the satellite at all times 
for a given orbital inclination, altitude, etc. This approach requires 
that every combination of orbital inclinations, sweep width, altitude, 
etc., must be investigated to have a complete data bank for estimating 
the frequency of observation or the time before an event is observed. 

This, by nature of the model's complexity, requires the use of a high 
speed computer and much computation time. For example, if each of three 
parameters is to be investigated for 10 different values, then 1000 separ- 
ate replications would be required. 

The author believes that an alternative approach, by using a simpli- 
fied analytical model, may prove to be just as useful for most appraisals 
of satellite surveillence systems. Appropriate analytical expressions, 
even if approximations to the actual situation, usually help to point out 
the important parameters, variables, etc. As a consequence, a simplified 
model like that discussed in this paper may point out some possible con- 
clusions not immediately apparent in a computerized model. For example, 
it was mentioned that the United States launched seven surveillence satel- 
lites during 1965 which were characterized by polar orbits. It will be 
obvious from the model to be presented in this paper that for reconnaissance 
purposes a polar orbit is extremely inefficient unless the polar regions 
themselves are of particular interest. 

Basic to the approach is the development of a fairly simple analy- 
tical expression for the likelihood of seeing any event of interest. It 
would be convenient if we could find the probability that the satellite 
passes over a chosen ground position during any randomly selected orbit — 



7 



but since the satellite’s position is completely deterministic we can not 
speak of the ’’probability” that the satellite passes over such a position. 
(For a selected orbit either the satellite passes over the position or it 
doesn’t.) However it does make sense to speak of the ’’fraction of orbits” 
in which the satellite passes over a selected ground position. This dis- 
tinction between ’’fraction” and ’’probability” is a fine one, but an impor- 
tant one. For example, if a satellite sees London every seventh orbit then 
it sees London on one-seventh of the orbits, but the probability it sees 
London on any randomly selected orbit is not one-seventh. 

The likelihood function to be derived will, then, describe the frac- 
tion of orbits that the satellite spends in ’’positive” orbits, where a 
positive orbit is defined to be one such that the satellite passes over 
the ground position in question on that orbit. If, for example, the ground 
position is Beirut, then the fraction of positive orbits would be the ratio 
of the number of orbits in which Beirut is within the sensors sweep width 
to the total number of orbits (per some time interval, e.g. the life of 
the satellite). If P represents this fraction, then, using reasonable 
assumptions, P can be shown to be a function of the satellite’s orbital 
inclination, the sweep width of the sensor, and the latitude of the point 
of interest. P is also a function of the satellite’s altitude, in that 
the sweep width is usually a function of the satellite altitude. 
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2. Formulation of the Likelihood Model 



The formulation of the likelihood model is facilitated by the 
following sketch. 




Satellite motion in general is confined to being great circle mo- 
tion. That is, the ground track of the satellite describes a great 
circle on a stationary globe. Controllable parameters are altitude, 
orbital inclination, and sweep width. 

For surveillence purposes, upper and lower bounds exist on satel- 
lite altitude. The lower bound is caused by the earth's atmosphere. A 
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satellite must be outside the atmosphere to remain in orbit. The upper 
limit is dependent on the resolution capabilities of the satellite sen- 
sors. The lower and upper bounds are presently about 100 and 300 miles 
respectively. The 300 mile maximum altitude means that the time for one 
orbit will be less than 1.4 hours. That is, while the earth rotates one 
revolution the satellite has orbitted about 17 times. This formulation 
will assume that altitude and sweep width are directly related and that 
knowing sweep width is equivalent to knowing satellite altitude. There- 
fore no further mention of altitude need be made. 

The orbital inclination angle -0- is equal to the extreme latitude 
attained by the satellite. Once launched, we assume & remains unchanged 
for the life of the satellite. 

The sweep width 2r is the width of the path on the earth which the 
satellite sensor sweeps out as it proceeds on its orbit. For photographic 
sensors 2r wuld be the width of the area photographed on one exposure, 
for radar it would be the length of one scan, scans being perpendicular 
to the direction of the satellite motion. 

We let o< represent the latitude at which an event occurs, and since 
we have no control over the location (or time) of events of interest, o< 
is an uncontrolled parameter. The ground position at latitude c< which 
is of interest will be designated as A, where A is either selected ran- 
domly or non-randomly , depending on the context. 

General Assumptions of the Model 

The following general assumptions are basic to the likelihood model: 

(i) The satellite under consideration remains in orbit for long 
periods of time, upwards of 1000 orbits. The life of a sat- 
ellite will be denoted by L. 
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(ii) An event of interest is equally likely to occur at A at any 
time. Given an event occurs during the time L, then we 
assume the random variable denoting the time of occurrence 
is uniform on (0,L). 

(iii) 9 cannot be varied after launch. 

(iv) The orbit is circular. 

(v) Every point on earth which lies within the sweep width of the 
satellite's sensor as the satellite passes overhead is seen 
with probability 1. 

(vi) The earth is a perfect sphere. 

(vii) The sensor points normal to the earth's surface at all times, 
and can see a distance r to each side of the satellite’s 
ground track. This varies with the actual sensor aspect in 
designed systems (in general the sensors are aimed to the 
side) but simplifies the trigonometry, without changing any 
basic concepts of the system. 

(viii) Sensors are either photographic or radar. In either case, 
the ensuing model assumes perfect visibility. (Brief men- 
tion is made in Section 4 of how one might incorporate the 
weather problem, and consequently the problem of general 
atmospheric attenuation, into the model) 

The Question of Accuracy 

The formulation to be presented does not yield extreme accuracy of 
results. Extreme accuracy seems unnecessary in the problem of initial 
system evaluation. For example, if we determine, through the model, that 
the mean time to next sighting of a ground position A is 2 hours when in 
actual fact it is 2.078 hours, then this probably would be of small con- 
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sequence in any decision made regarding the system. Certainly it would 
be misleading to present results implying accurately to six or seven fig- 
ures when basic assumptions preclude accuracy to more than two or three 
figures. Computer formulations specifically designed for the surveillence 
evaluation problem may claim accuracy to six figures. [3} These figures 
may or may not be that accurate when compared to an actual satellite in 
orbit, but any initial evaluation of a satellite system would probably have 
remained the same had the figures been rounded to two significant figures, 
which is approximately the accuracy inherent in the proposed model. 



We will denote the ’’fraction of positive orbits” by P(# , o< , r) . 
However P(& , o< , r) will often be referred to as the likelihood function 
P. As the notation implies, P depends on three parameters: orbital in- 
clination 0, latitude o< of the ground position A to be observed, and the 
sweep width 2r. 

The analytical expressions for P(£ , , r) are derived in Appendix 

I and are given below. The derivation is based on the assumption that P 
is equivalent to the ratio of the sweep width measured along a latitude 
to one-half the circumference of the latitude. (i.e. TS in the previous 



The Likelihood Function 



sketch) 
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where S = r and R is the radius of the earth. 

R 

Figure 1 illustrates the behavior of P as a function of <x for 
several values of 6 when r = 150 nautical miles. The figure shows 
that P attains its maximum when o( = . This fact can be con- 

firmed by an investigation of (!)• 
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Figure 1. P(0,«,r) as a Function of o< for Fixed Values of r and 0 



3. Applications- 



Mean of Time Until Next Sighting 

One of the present parameters used for appraising a satellite sur- 
veillence system is the mean time between sightings.* From our model the 
mean time between sightings is simply -Is_ . However if location A has 
suddenly become of interest we probably would not care about the time 
lapse between the last sighting of A and the next one, but rather on the 
time lapse from the present until the next sighting. These tw parameters 
are by no means equivalent. 

In this section we will develope several estimates for the mean time 
from a given instant until the next sighting of a specified point A. In 
these developments, we will assume that phenomena of interest are equally 
likely to occur at A at any time during the life of the satellite. 

If T is the time lag between occurrence of an event at location A 

Tj 

and the first sighting of A after the occurrence, then we are looking for 
the mean value of T L , denoted by M(T t ) . An exact expression for M(T L ) is 
difficult to obtain, but various estimates can be made. 

To see why the exact expression is difficult to obtain analytically 

consider the following sketch. 




^■For an example, see ClJ • 
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If we consider an upward (downward) pass through the latitude c< to be one 

such that the satellite is moving from o< to a more northerly (southerly) 

latitude, then the satellite can see point A on any given orbit in one of 

two ways. Either it sees A on an upward pass or on a downward pass. The 

time between an upward and downward pass depends on the latitude o< . If 

T* deontes this time and T is the time for a satellite to make one com- 

o 

plete orbit, then the discrete points in time at which sightings can pos- 
sibly be made occur in a sequence T^, T2, T^, where the times 

between sightings are T* , T Q - T* , T* , T q - T* , . . . . 

In reality, sightings do not occur in a neat sequence T^, T2,..., T^, 
... as described above. That sequence merely tells when sightings are pos- 
sible . The actual sightings may occur in sequences which are exceedingly 
difficult to describe analytically. For example, as indicated by the 
previous sketch, a sighting occurs when on an upward pass the satellite 
passes over point P while A lies between Q and Q*. Since the distance 
QQ t is on the order of several hundred miles, the speed of the earth 1 s 
rotation makes it possible for A to transit QQ* in much less time than it 
takes the satellite to make a complete orbit. Thus the satellite may have 
just crossed o< while A was just west of Q, and A will pass Q* before the 
satellite makes another upward pass. A "sighting possibility" would thus 
be missed. 

Rather than attempting the difficult task of predicting the actual 
sequence of positive orbits, we will develope estimates for M(T^) . If 
tests show these estimates to be reasonably accurate, then they should be 
almost as useful to the decision maker as the true value. Let us consider 
three possible estimates. 

Estimate I; If we assume actual sightings occur at equally spaced inter- 
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vals of time, then the spacing can be obtained from the likelihood function 
P. The assumption that events of interest occur with equal likelihood at 
any time during the life of the satellite implies that if an event occurs 
during the time (T^, T^ + p then it is equally likely to occur at any in- 
stant in that interval; that is, the conditional time of occurrence T is 

T 

uniform on (T^, T^ + ^). Since the length of this interval is , M(T^) , 

an estimate for M(T^) , is 



M(T l ) = 



2P 



Intuition suggests that this may be a good estimate near the equator 
where the time between an upward and downward pass, and therefore the time 
between possible sightings, is the constant %T q . Since the times of pos- 
sible sightings are equally spaced, the time between actual sightings may 
tend to be equally spaced also. But it seems equally intiutive that as A 
moves away from the equator the estimate will become less reasonable. 
Estimate 2 : Arbitrarily define a probability model for the orbiting sat- 

ellite, wherein P, instead of being the fraction of positive orbits, is 
defined to be the probability of a positive orbit. Further assume that 
all orbits are independent of each other. This model is then equivalent 
to one with Bernoulli trials with probability of success equal to P. Fur- 
thermore the number of trials to first success for such Bernoulli trials 
is a geometric random variable with mean — i— ; where P is the probability of 
success on each trial. Equivalently, in terms of our model, the nunber of 

orbits to the first success (positive orbit) is - 1 — and therefore the 

T 

estimate for the time to first positive orbit would be . 

Intuition suggests that this estimate should be good when o( ap- 
proaches 0, but not so good near the equator. The time between an upward 
and downward pass approaches zero as Ot-*'©*, and between downward and up- 
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ward approaches T . Thus the time till possible sighting will be twice as 
long as that at the equator. 

Estimate 3 : Since neither of the above estimates seems reasonable for all 

values of o< , we return to the physics involved to see if an assumption 
can be made which provides aubetter estimate. 

If a ground position A is at latitude c* we have said it can be seen 
during that orbit at one of two times, it can be seen during the upward 
pass or the downward pass. It is relatively easy to obtain an exact ex- 
pression for T* , the time between the upward and downward passes. Knowing 
T q and T* we can describe the sequence of times at which sightings are 
possible. We also know that the fraction of positive orbits is P. 

Using these two pieces of information we assume that actual sightings 
occur in a sequence such that the time between sightings is proportional 
to the time between possible sightings, where the constant of proportion- 
ality is such that the fraction of positive orbits is P. We get (see 
Appendix II) 

M (r L ) - cof & fan <x)] 2 -h 

£tt - C05~ / ( : cof ■& fan c< )J Z J ( 3 ) 

A T 

Notice that for 0( = 0,M(T t ) becomes — _ which is identical to 

L 2P 

a T 

estimate 1, and for o( = 0 , M(T 7 ) becomes ° which is identical to 

L P 

estimate 2. Thus estimate 3 varies with o( and also fits what intuition 
suggests should be the case at the extremes of the range of o( . Further 
evidence of this estimate’s reasonability is given by Appendix III where 
M(T^) is compared with results obtained from a computerized simulation 
developed by the Naval Missile Center. £3] 
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Location and Number of Readout 



Stations Required for a Specified Level of Contact 
Any study of the effectiveness of a satellite surveillence system 
must include an appraisal of problems of transferring ’the information 
from satellite to evaluation centers located on earth. To consider a sat- 
ellite which orbits for long periods of time we should discuss readout 
stations. Readout stations are radio receiving stations capable of re- 
ceiving signals from a satellite whenever the satellite is within a cer- 
tain distance of the readout station. We will refer to the circular area 
wherein such reception is possible as the transmission area . As the sat- 
ellite passes through a transmissidn area it transmits information accumu- 
lated since the last pass through a transmission area. The author knows 
of no present study which analyzes the optimum location of such stations 
or the number needed to ensure adequate transmission time, therefore it 
seems worthwhile to consider this question. 

P(£, o< , r) can be used to determine the percentage of time a satel- 
lite will be in contact with a particular readout station. The procedure 
we will use is: 

1. Find the diameter of the transmission circle within which 
radio signals from the satellite are receivable. 

2. Estimate M(t) , the mean time a satellite takes to cross the 
transmission circle given it intersects the circle. 

3. Determine the value P(0*, oC , r) , where o( is the latitude 
of the transmission circle and r is the radius of the transmission circle. 

4. Estimate the mean time per orbit available for transmitting, 




Knowing the rate of transmission of information, a reasonable esti- 
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mate for the number of readout stations needed to maintain a specified 
level of information transmission can be obtained. 



To find the diameter of the transmission circle consider the fol- 




In this sketch h is the height of the satellite above the earth’s sur- 

/ 

face, and a is the minimum reception angle above the horizon. 

From the law of sines, 

Sin (an<j!& A&C) _ 3 j n (ai-g)_ _ cos a 

R " R R+H 

/ here -fore > Anjk ABC = s/n ' 

W, cos' Y - a. 

Thus, the diameter of the transmission circle, d, is 

J= * . W 
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Thus we see that a readout station can receive signals from the satel- 
lite if the satellite is within a circle with diameter 2(R + h)0 with 
center at the readout station. If we no# let r in the expression for 
P(*0> > r) equal 1(R + h)0 , then we will have the expression for the 

fraction of orbits wherein the satellite passes within the transmission 
circle. 

To find the mean time available for transmitting given the satellite 
passes through a transmission area, we will assume that the distribution 
of paths which the satellite makes as it crosses the transmission area is 
uniform; that is, the satellite is equally likely to cross the circle on 
any chord. Finding the mean length of path based on this assumption is not 
difficult. If the radius of the transmission circle is r then the mean 
distance travelled while in the circle is approximately . (See Appen- 

dix IV) . 

From the above results, we can obtain an estimate for the mean time 
per orbit spent in the transmission area of a single readout station. It 
is 

M(t) = P(6, * , r); r = (R + H) 

V = ground velocity of satellite. 

Because P is maximized when o( = ■©* - £ and M(t) is linear in P, we 
can maximize the mean transmission time per orbit by placing the readout 
station at latitude -0* - & . 

We next consider the question of how many readout stations are re- 
quired to attain at least some specified mean transmission time per orbit. 
Wa can reasonably assume that all stations would be placed close to lat- 
itude 4-8 . Wfe also assume the stations are spread out on that lati- 
tude to avoid a satellite passing through two stations in rapid succession 
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and then not passing through another for many orbits. Under these assump- 
tions a reasonable estimate for the number of readout stations N required 
to average t 1 minutes of transmission per orbit is 



N = 



t 

t 



If a value of t’ is designated as the minimum mean time per orbit needed 
in order to receive some specified fraction of the total information ac- 
quired by the satellite, then N tells us how many readout stations are 
required. Such information wuld be needed to estimate the cost of a pro- 
posed satellite system. 

Estimating the Number of Intersections of a Transitting Ship 

by a Satellite 

To find an estimate for the number of intersections of a moving point 
on earth, consider a ship going from point a to destination b, where a 
and b are at latitudes A and p respectively. Let us assume: 

1. The ship proceeds at constant speed s. 

2. The ship moves on great circle path. 

3. The distance d from a to b on this path is known. 

4. The track of the ship does not leave the area of satellite 

coverage. 

First, to avoid complicated expressions derived from differential equa- 
tions giving the instantaneous latitude of the ship as it proceeds, we 
elect to use a constant, % , to represent the time averaged mean latitude 
of the ship; that is , ft is an estimate for the latitude which, if used to 
replace o< in P(&, o( , r) , would give the same result as if we continu- 
ously evaluated P at each instant of time during the ship's transit. IS is 
defined as . 

Y A +<3 

O - 2 
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Since the ship's track is a great circle track and the ship's speed is 
constant, this is a lower bound for the ship's time averaged mean latitude. 
A conservative estimate for the number of intersections will consequently 
be obtained. 

The derivation of the estimate of the total number of sightings for 
a ship transitting from a to b under the above assumptions is the follow- 
ing: 

Let h = satellite altitude 

t = ship’s transit time from a to b 
v = satellite's speed 
g = gravitational constant. 



Then 



= / Z I —— , and t = — fL 

y (R + h) s 



Now 



T = 2 TT (R + h) 

O v 



and thus N, the number of orbits during the transit is 

N = 4 — . 



If I is the number of intersections during the ship's transit then an 
estimate for I is 

A _ 

I = NPC-0-, 2T , r). 



Example: If the ship goes 1000 miles from 10° to 20° at 12 knots, the 

average number of sightings by a 200 mile high satellite with orbit 
9 = 50° and r = 150 will be 

I = P • N = NP(£,tf , r) = (63. 2) ( .035) = 2.2. 
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4. Areas for Further Study. 

Having derived an expression for the fraction of positive orbits one 
can investigate various areas other than those already mentioned. The 
following are areas which seem to be of obvious interest in the evaluation 
task. 

1. Finding the optimum orbit based on some measure of effectiveness 
is a problem of considerable interest. Figure 2 shows graphs indicating 
the distribution of major nautical activities as a function of latitude. 
Further work could be accomplished correlating P and these distributions 
to optimize the orbit. For example, if the measure of effectiveness is 
"number of contacts per unit time" then we could proceed as follows. Des- 
cribe a graph of figure 5 by a function f(oO, where f is either a discrete 
or continuous function approximating the graph. If we consider southern 
latitudes to be negative, then we can let 

* 

F(o() = f ffa'Jc/*' or F(<x) - Z2 {(*') ) 

depending on which is applicable. Then, let 

G (& L > r) = j' P(&><*,r) c/Ffo() , f/tfj 

If G is the maximum of G with respect to 0^, then that 9^ which yields 
G is the optimal orbital inclination for our chosen measure of effective- 
ness. 

The distribution of naval activities with longitude may be combined 
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Figure 2. Distribution of Various Naval and Maritime Activities 
with Respect to Latitude. [2j 
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with the distribution for latitude. The optimum orbit for observing speci- 
fied regions may then be quite different than if latitude distribution 
only are considered. 

2. In studying the effectiveness of a satellite as a weapon delivery 
vehicle, we need only replace the half sweep width r by the lethal range 
r f of the weapon carried by the satellite. Then, if the ground position 
of the satellite is within the lethal range r f of the target, the weapon 

A 

could be dropped. P and M(T^) would have the same form as above, but 
MCT^) would be the mean time until a weapon could be dropped on a target 
after the decision to destroy the target had been made. One might cor- 
relate this with the possibility of inflight satellite course changes. 

3. Investigating countermeasure techniques may prove of interest. In 

this respect figure 1 gives some indication of P and its sensitivity to 
latitude in the extreme regions near 0. Knowing the inclination of an 
enemy’s orbiting satellite and correlating that information with the impli- 
cation of figure 1 would indicate ship routing procedures we could use to 
avoid detection. We obviously would try to route our ships to stay clear 
of the latitudes near $ - £ , where P is largest. 

4. When data is available on global cloud cover distributions, P can 
be modified by the probability of cloud cover and/or darkness obscuring 
the sensor’s vision (as in the case of photographic sensors). As a con- 
sequence, P would become closer to being a probability than a fraction. 

As a result of cloud cover being a function of longitude as well as lat- 
itude, P would also have to be dependent on longitude. The problem of 
predicting cloud cover is obviously complex and any model developed to 
describe the distribution of cloud cover would probably have to be a com- 
puter simulation. 
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APPENDIX I 



DERIVATION OF THE FRACTION OF POSITIVE ORBITS 
The function , r) is defined to be the fraction of 

positive orbits to total orbits given a satellite in orbital inclina- 
tion 0, a sensor sweep width 2r, and a ground position at latitude c<. . 
The basic rationale for finding P is to consider the circle on the 
sphere made by latitude ci , and to find the ratio of the path width of 
the satellite measured along latitude u to one-half the total circle. 
Then if the ground position A is on latitude o£ , the fraction of times 
A is seen will be proportional to this ratio. 

To obtain an analytic expression for P based on this rationale 
consider the following sketch, in which capital letters denote points 
on the sphere, lower case letters denote angles of spherical triangles, 
and Greek letters denote sides of spherical triangles, a side being the 
usual radian measure of the central angle subtended by that side. 
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If we could determine the length of BB*, then divide BB 1 by FG, (one 
half the circumference of the circle made by latitude ) we would have P. 
However it is not necessary to find BB 1 . We get the required results by 
finding, instead, angles a and c, subtracting these from tt to get angle 
b, then dividing b by TT . 

To find angle a we will make use of the fact that triangles AND and 

CBD are spherical right triangles* Furthermore, side Y* of AND is 0. 

Thus if we can find side A then Napier 1 s Rules for spherical right tri- 
angles will give us angle a. A similar approach is used to find angle c. 

We first review Napier's Rules. Given a spherical right triangle as 
shown in the sketch below, where <X , (3 , and represent the central 
angles subtended by the respective sides of the triangle, then the rules 
can be summarized as follows: 

Take the five parts, excluding the right angle, and con- 
sider them to be in the circular arrangement shown. If 
"co- 11 means "complement of" then: 

(1) the sine of the middle part equals the product of 
the tangents of the adjacent parts; 

(2) the sine of the middle part equals the product of 
the cosine of the opposite parts. 





A 
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To find angle a we first note that g is a right angle since A is 
the northernmost point of the orbit, and h is a right angle by con- 
struction. The sweep width is 2r and the central angle subtended by r 
la then 




From Napier’s Rules, 

sin i ~ cos(co-S) cos(co-e) 



and therefore 



sin § 



sin 
stn e. 



Similarly, 



But 



sin X 



sin Y 3 
sin e 





(/z) 

(/ 3 ) 



From ( 11 ), (12), and (13), 

sin Y 5 sin$ 



s in X - 



sin 'f cos (v.+ a) 



stn 



sin 8 



. 0 +) 



After some manipulation of (14), the expression for X is 




Os) 



Again using Napier’s Rule, 

COS CL - Al/7 'f coi X . (i6) 



From ((5) and (16) , 



a = 



-[sec ?[■ 



sin £ + cos(£^)_sin_ )° 77 /,y\ 

s tn(£ + k) JJ- [ ' 
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To find angle c consider the following sketch. 



N 




Due to symmetry between northern and southern hemispheres, B*G equals 
B* *G* ' ; therefore angle c equals angle k. But angle k obviously equals 
angle f; therefore to find c we need only find f . Now f can be found in 
the same way we found angle a by merely using spherical right triangles 
AND' ’ and C'^'D* * instead of triangles AND and CBD. Therefore we can 
use equation (17) directly by replacing ( \ + £ ) with ( V + §')• But 

( A ’ ♦ £’) must equal 7T -( A + § ) because (§* + A’) equals ( 7T - l) ) 

and equals ( X + § ) by symmetry. Therefore, after the appropriate 
changes in (17) , we get 



C = cos' ■ L sec 



r'fsec t[- 



sin £ + cos ( tt— % -X) sin Y* 



sin (ir-g-X) 



-]]■ 



( 13 ) 



\ 
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We are now ready to obtain an expression for PC’d' , <* , r). We 
said P would be equal to , where b is obtained from 



b = 7T- (a + c) . 



(/V 



Using the relationships 

§■»* = f- « , 

v = ¥-*> 

■and P (&>**•> r ) - 



in equations (17), (18), and (19), we obtain / 



pc *. - / - ¥[^"[0 ' ( sm *+ si ” * + 

cos-'[ f^( s/nS ~ S '”‘ < cos *)]}■ ( Z °) 



On a reexamination of the derivation we realize that (20) is valid only 
for £■ - $ . Because angle a is zero for <* > & - £ , (20) must 

be modified by dropping the first term in for 

Finally P must be zero for <s( > -fr + £ • Therefore the expression for 
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P over the range O £ <* - is 

Cm 



PiW)=< 



/-fr[ --°s'[^r( s/n s *™"* cos ^ 

[S^-fsm S - cos <§]-, &i J s 

COS lcc>i<n[ -rj / 0£ o( < -O—S 



o ; *+s *«*■£• 



L 



The behavior of P over the range 0 is illustrated in 

figure 1. 
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APPENDIX II 



Derivation of Estimate 3 for Msan Time to Next Sighting 
Consider the following sketch: 

/V 




Wa will consider our line of sight to be always across the plane 
of satellite rotation, and let the earth turn relative to that plane. 

Let D = time for single revolution of earth in hours 

D = time in hours for earth to turn far enough so that A (fixed 
on earth) moves twice the distance QR. (Thus D' is the 
time between a possible sighting on an upward pass, and a 
possible sighting on a downward pass.) 

T q = Orbital period in hours 

From the sketch, 

D’ = - arc <? R D = an gle g D 
arc PR 7T 

But, by Napier's Rule, 
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cos g = cot ^ tan o< . 



Therefore D’= cos" 1 (cot » tan « ) ^ 

rr 

We want to find an estimate for the mean time to next sighting of 
location A, given A has suddenly become of interest. We are unable, with- 
out great difficulty, to describe the sequence of actual sighting times. 
Therefore we would like to define a sequence of times which we can use to 
approximate the sequence of actual sightings, and upon which we can base a 
reasonable estimate for M(T L ) . , *» 

We will assume that the sequence selected should satisfy two condi- 
tions. 

Condition 1 . The sequence must satisfy the condition that the 

mean number of orbits between positive orbits in n, where n = — L- . 

P 

Condition 2 . The sequence should be proportional to the sequence 
of possible sighting times. 

It is obvious why condition 1 is desired. Condition 2 is specified be- 
cause it will require the sequence of assumed actual sightings to be 
spaced in a way which suits our intuition. Thus, if the possible sight- 
ings occur in a sequence such that the time between possible sightings 
has the sequence 

1 hr, 23 hr, 1 hr, 23 hr, 

then we would intuitively expect actual sightings to be spaced such that 
usually two sightings would occur relatively close together, followed by 
a long period with no sighting, then two more sightings relatively close 
together, etc. 

Therefore if the times between possible sightings has the sequence 
Z = D’ , D-D f , D* , D-D 1 
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then we want a constant K such that the sequence KZ satisfies our two con- 
ditions. There is obviously only one sequence which satisfies these two 
conditions simultaneously . * Therefore our next task is to find that se- 
quence. 

To find the sequence KZ we need only find a K which, when multiply- 
ing Z, satisfies condition 1, because any sequence KZ where K is constant 
will automatically satisfy condition 2. Now the sequence Z ’’cycles'* on 
every second sighting possibility with a cycle time of D and therefore 
there are two possible sightings in a cycle of Z. Wfe would like KZ to 
’’cycle” on every second sighting also, but the cycle time must be equal to 

2nT Q , since there are to be 2n orbits per every two sightings. Therefore 

2nT 

to get KZ we must multiply Z by -2— . Thus: 

/r= ^2— • 



Wa then have 



KZ= [ 



» ' ~D~ 



(P-D‘), ZjlHd: izUfo-e') ■■■]. 

D D > J 



then, 



To show that this sequence satisfies condition 1 we consider a single 

O 



cycle n 7 q r,' The mean time between sightings is, 



(p-p'J] = nT 0 



and therefore the mean number of orbits between positive orbits is n. 



1 We have fixed the mean time between points in the sequence and 
the "cycle" of the sequence. 
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Having defined KZ, the sequence of times between assumed actual 



'ic 

sightings, we can say that assumed positive orbits occur at times T^, 
* * 

T 3 , T 3 , where 




f n* * i odd 

' r* y 2n Z(£z£!) i i e™ . 



/•e. 



(F, 5 



x -r * 

> '3 



T* 

'v > 



"} 



[ZnZ (g), 2n T a , 2«T 0 + 2 ”(£'), -] 



To estimate M(T^) on this basis consider the following diagram: 



-+-T . 1 


-2nT 0 

| n * 7V * 


. T — .1 


1 




'6, \ 








where Tj indicates the time of the jth assumed sighting and the time be- 
tween and T^ + 2 *- s defined to be one cycle. Then assuming the time of 
occurrence of an event of interest at A is uniform over the satellite’s 
life, we can say that the event of interest is equally likely to occur 
at any time T contained in (T^, T^ +3 ) • 

We are trying to estimate the mean, time to next sighting given an 
event of interest has occurred at A. However, since each cycle is iden- . 
tical, we can equivalently find the mean time to next sighting given 
T € (Tp • The density of T, where T is the random variable denoting 
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the time of occurrence of the event of interest given T € (T? , T? +2 ) , 
is then 



f r (t) 



ZnT 0 



Therefore 



a f , Tq. fmean h'me irom T /o 7^*/ 7 

M(rJr ^Ls^ recr.rj)] 



+ 




Substitution of the appropriate expressions for the 



Tj's in (22). results 



in 




and, upon substitution for D' from (21), 

A(T l ) = -J^J[c0S- '(ce+e-hn*}! 2 V- 

- Cc ?s~‘ (ce+& -t-an I x)Jj. (23) 



Thus, we obtain a general estimate for M(T^) as a function of ■& , u , and 
r. 



\ 
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It is of interest to note that when o( = 0, (23) reduces to 



M(T l ) 



12 - 

2P ’ 



and when o( =&, (23) reduces to 



M(T l ) 




These limiting forms are identical to estimates 1 and 2 respectively, 
which is a desirable result, as estimates 1 and 2 seem appealing at these 
limits. 
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APPENDIX III 



Comparison of Theoretical and Numerical Estimates for M(Tl) 

Our analytical results are compared with simulated results for the 
following data: 0 = 50°; r = 150n.nu; T Q = 1.6 hrs; 105 orbits. 



cx 



0 

20 

40 

50 



Analytical Results 
32 hrs 
22.7 hrs 

16.5 hrs 

12.5 hrs 



Simulated Results 
35 hrs 

22.5 hrs 
16.7 hrs 

10.5 hrs 



The simulated results were obtained from a Stanford Research Insti- 
tute study which computed the times of intersection for each of some 1800 
fixed locations. [2] To obtain M(T^) from these times of intersection 
for each location the occurrence of an event was assumed to be uniformly 
distributed over time during the 105 orbits considered. M(Tj) for loca- 
tion j was then found by the formula 



n j 






where o^* \ tlj are the times of intersection, jdj is the total 

time considered for location j(/^j = time between the first and last in- 
tersections for j) . This formula is easily verified as a formula for the 
time to sighting of an event which occurs uniformly in (0, /^j) when times 
of sighting occur at times cX, , <X 2 } , Then the overall value of 

M(T t ) was found by averaging the M(T-). Since locations were scattered 
L J 

throughout various latitudes only those within a small band were considered 
to find M(T l ) for that band. For example, all locations between 39° and 
41° were used to find M(T ) for 40° latitude. In the case where two cities 
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were so close together that considering both would create bias, (Los 
Angeles and Long Beach) only one was used. 

The results are extremely close with the notable exception of the 
range where we would expect our estimate to yield the best results, i.e. 
at - 0 = 50°. The difference here is most likely due to the fact 
that the model in this paper assumes the satellite sensor "looks straight 
down" whereas the computer model has the sensor's field of view skewed 
toward the north side of the path. This difference between our model 
and the computer model can be seen from the following sketch. 




The satellite in the latter case, when it is over latitude 50°, 
could have its sensor centered on a latitude as high as 54°. While a 
four degree error is negligible at moderate latitudes, as approaches 
0 the difference becomes quite large. If the fraction of positive orbits 
from the computer simulation is computed and used in our formula instead 
of P(9, o< , r) , then mCT-jO for o< = 50° is found to be 10.7 hrs. There- 
fore the variation between our results and those of the computer simulation 
appears, to be caused by the difference in fraction of positive orbits due 
to slightly different physical models. 
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APPENDIX IV 



Derivation of Mean Distance to Cross Transmission Area 

We want to find the mean distance X that a satellite travels when 
it transits a transmission area of a ground readout station. We can 
estimate this distance by assuming the transmission area is a plane cir- 
cle, and that the track of the satellite is a chord through the circle 
with the distance from the chord to the center being uniformly distribu- 
ted on (0,r). To facilitate the discussion, consider the following • 
sketch. 



In terms of the sketch we want to find E(X), the mean of X, when 
the random variable Y is uniform on AC, that is, 




X 




From the sketch we get 
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and, because Y is uniformly distributed, the random variable X* has 
the distribution function: 



F x .(*)= /-Pr(V*f^f) 

= /- o*X*r. m 



But X = 2x' ; therefore the distribution function of the random variable 
X is 




OiFiZr: (ZS) 



• It is then easily shown that 



£(X) = 





Equation (26) was obtained assuming a plane circle, whereas the actual 
transmission area is a segment of a sphere. The' result is a reasonable 
approximation, however, because the transmission areas are small enough 
that the curvature of the earth is negligible. 
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